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Abstract. We investigate the analytic structure of solutions of 
non-relativistic Schrodinger equations describing Coulombic many- 
particle systems. We prove the following: Let ip(x) with x = 
(afi, . . . ,Xff) S R 3N denote an TV-electron wavefunction of such 
a system with one nucleus fixed at the origin. Then in a neigh- 
bourhood of a coalescence point, for which x\ — and the other 
electron coordinates do not coincide, and differ from 0, ip can be 
represented locally as = ip^ijx.) + \x\\ip( 2 \-x.) with ip^, ip( 2 ' 
real analytic. A similar representation holds near two-electron coa- 
lescence points. The Kustaanheimo-Stiefel transform and analytic 
hypocllipticity play an essential role in the proof. 



1. Introduction and results 

1.1. Introduction. In quantum chemistry and atomic and molecular 
physics, the regularity properties of the Coulombic wavefunctions ip, 
and of their corresponding one-electron densities p, are of great impor- 
tance. These regularity properties determine the convergence proper- 
ties of various (numerical) approximation schemes (see [26, 2, 3, 31, 
32, 33] for some recent works). They are also of intrinsic mathematical 
interest . 

The pioneering work is due to Kato [20] , who proved that tp is Lips- 
chitz continuous, i.e., if) G C 0,1 , near two-particle coalescence points. 

In a series of recent papers the present authors have studied these 
properties in detail. In [7] we deduced an optimal representation of 
ip of the form if) = JF$ with an explicit T G C 0,1 , such that $ G 
C 1 ' 1 , characterizing the singularities of ip up to second derivatives; see 
[7, Theorem 1.1] for a precise statement. In particular, T contains 
logarithmic terms which stem from the singularities of the potential 
at three-particle coalescence points. This characterization has been 
applied in [8] and [9] in the study of the electron density p and (in the 
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atomic case) its spherical average p close to the nuclei. Real analyticity 
of p away from the nuclei was proved in [6] ; see also [4] , [5] . 

In this paper we derive a different representation of ip which com- 
pletely settles its analytic structure in the neighbourhood of two-par- 
ticle coalescence points. The Kustaanheimo-Stiefel transform (KS- 
transform for short) and analytic hypoanalyticity of a certain degener- 
ate elliptic operator are crucial for the proof. 

We start with the one-particle case. 

Theorem 1.1. Let Q G M 3 be a neighbourhood of the origin and as- 
sume that W^ 2 \ F^\ and are real analytic functions in f2. 
Let 

H = -A + ^-^ + W^, (1.1) 

\x\ 

and assume that ip G W l,2 {VL) satisfies 

pi 1 ) 

H<p=— + F& (1.2) 
\x\ 

in Q in the distributional sense. 

Then there exists a neighbourhood Q G Q of the origin, and real 
analytic functions : VI — > C such that 

cp(x) = <pW(x) + \x\v {2 \x) , xett. (1.3) 

Remark 1.2. Theorem 1.1 is a generalization of an almost 25 years old 
result by Hill [16]. The present investigations were partly motivated 
by this work. Hill considered solutions to 

( - A - ||r + V {1 \x) + \x\V {2 \x)) V = , (1.4) 

with and real analytic near the origin, and proved that (p satis- 
fies (1.3). The statement (1.3) can easily be seen to hold for Hydrogenic 
eigenfunctions. These have the form e~ l3 ^P(x) for some (3 > 0, where 
P{x) can be written as linear combinations of polynomials in \x\ times 
homogeneous harmonic polynomials. In particular, Hill's result implies 
that ip satisfies (1.3) near the origin for a one-electron molecule with 
fixed nuclei, one of them at the origin. 

Remark 1.3. Hill's proof is rather involved. Our proof is quite differ- 
ent, also not easy, but has the advantage that it can be generalized to 
treat the Coulombic many-particle case; see Theorem 1.4 and its proof 
below, and also Remark 1.6. 

The proof of Theorem 1.1 uses the KS-transform (see Section 2 for 
the definition). This transform was introduced in the 1960's [25] to 
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regularize the Kepler problem in classical mechanics (see also [24, 30, 
22]) and has found applications in problems related to the Coulomb 
potential in classical mechanics and quantum mechanics, see [1, 10, 13, 
14, 15, 19]. The KS-transform is a homogeneous extension of the Hopf 
map (also called the Hopf fibration), the first example of a map from 
§ 3 to S 2 which is not null- nomotopic, discovered in the 1930's [17]. For 
more on the literature on the KS-transform, see [22, 14]. 

We move to the iV-particle problem. For the sake of simplicity we 
consider the atomic case and mention extensions in the remarks. Let 
H be the non-relativistic Schrodinger operator of an iV-electron atom 
with nuclear charge Z > in the fixed nucleus approximation, 

* = E(-A*-ifi)+ E ra =: - A + v '' (L5) 

j=l 1 31 l<i<j<N 11 n 

Here the Xj = (xj t i,Xj j2 , 2^,3) G M 3 , j — 1, . . . , N, denote the positions 
of the electrons, and the Aj are the associated Laplacians so that A = 
Y2f =1 Aj is the 3iV-dimensional Laplacian. Let x = (xi, x%, ■ ■ ■ , x^) G 
WL 3N and let V = (Vi,...,Vjy) denote the 3A^-dimensional gradient 
operator. The operator H is bounded from W 2,2 (M. 3N ) to L 2 (M. 3N ), and 
defines a bounded quadratic form on W 1,2 (R 3N ) [21]. We investigate 
local solutions ip of 

Hip = E^j , EeM, (1.6) 

in a neighbourhood of two-particle coalescence points. 
More precisely, let S denote the set of coalescence points, 

N 

S := jx= (x!,...,x N ) G R 3N j Y[\xj\ '[I \xi -Xj\ = oj . (1.7) 

3=1 l<i<j<N 

If, for some VL C ~R 3N , ip is a distributional solution to (1.6) in fl, 
then [18, Section 7.5, pp. 177-180] ip is real analytic away from E, that 
is, if) G C^iVt \ E). 

Let, for k,£ G {1, . . . , N},k ^ £, 



E'. 



{x G M 3JV j [] \ Xj \ J] \xi - Xj\ = 0} , (1.8) 



j=l,j^k l<i<j<N 
N 



\ t : = |x G M 3JV j |x* - x,| = 0} . (1.9) 



i=i i<i<j<N,{i,j}^{k,e} 
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Then we denote 

£ fc :=£\£' fc , £ M :=£\£' M (1.10) 

the two kinds of 'two-particle coalescence points'. 
The main result of this paper is the following. 

Theorem 1.4. Let H be the non-relativistic Hamiltonian of an atom, 
given by (1.5), let Q C M. 3N be an open set, and assume that ip G 
W 1,2 (fi) satisfies, for some E G K, 

Hit = Eib~ in Q, (1.11) 

in the distributional sense. Let the sets £ fc and lu k ^ be given by (1.10). 

Then, for all k G {1, . . . , N}, there exists a neighbourhood Q k C Q 
of fid £fe 7 and real analytic functions ij)^,^)^ : fife — * C such that 

V>(x) = tf } (x) + M^ 2) (x) , x G Vl k , (1.12) 

and for allk,£ G {1, . . . , iV} 7 k ^ t, there exists a neighbourhood fl k ,e C 
Q o/finSy, and rea/ analytic functions ip^},^] '■ Qk,e — > C suc/i £/ia£ 

V»(x) = ^( x ) + kfc - ^IV'mW > x e • (i- 13 ) 

Remark 1.5. The proof of Theorem 1.4 again uses the KS-transform. 
Due to the presence of the other electron coordinates we are confronted 
with additional problems. We have to deal with degenerate elliptic 
PDE's where the corresponding operators (of Grusin-type) are analytic 
hypoelliptic, see [12]. 

Remark 1.6. Theorem 1.4 extends in the obvious way to electronic 
eigenfunctions of Hamiltonians of iV-electron molecules with K nuclei 
fixed at positions (R±, . . . , Rk) G R 3 , given by 

N K 

" = £(-A J -£ R # 1? ) + £ J^-y (1-14) 

j=l l=\ 1 3 11 l<i<j<N 11 Jl 

Furthermore we can replace in (1.14), as in Theorem 1.1, the poten- 
tial terms by more general terms, and allow for inhomogeneities. 

For instance, the result holds for general Coulombic many-particle 
systems described by 

n A 

H = H-^t + Yl %(^-^), (i-is) 

j=l 3 l<«<i<n- 

where the rrij > denote the masses of the particles, and v i3 - = \xi — 
x j\~ 1 + v ij w ^ n v ij\ k = 1,2, real analytic. 
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Remark 1.7. In separate work we will present additional regularity 
results (not primarily for Coulomb problems) obtained partly using the 
techniques developed in the present paper. 



2. Proofs of the main theorems 

As mentioned in the introduction our proofs are based on the Ku- 
staanheimo-Stiefel (KS) transform. We will 'lift' the differential equa- 
tions to new coordinates using that transform. The solutions to the new 
equations will be real analytic functions. By projecting to the original 
coordinates we get the structure results Theorem 1.1 and Theorem 1.4. 

In the present section we will introduce the KS-transform and show 
how it allows to obtain Theorems 1.1 and 1.4. The more technical 
verifications of the properties of the KS-transform and its composition 
with real analytic functions needed for these proofs are left to Sections 3 
and 4. 

Define the KS-transform K : 1R 4 -> R 3 by 

(vl-yl-vl + vt\ 

K(y)= 2(y m -y m ) . (2.1) 

\ 2(ym + vm) j 

It is a simple computation to verify that 

\K{y)\ := \\K(y)\\ R3 = \\y\\ 2 R4 =: \y\ 2 for all yeR 4 . (2.2) 
Let / : IR 3 — > C be any C 2 -function, and define, with K as above, 

f K ■ K 4 - C , f K (y) := f(K(y)) . (2.3) 
Then for all t/GR 4 \ {0}, (see Lemma 3.1), 

(^f)(K(y)) = ^Af K (y). (2.4) 

2.1. Proof of Theorem 1.1. Assume (p G W 1,2 (Vl) satisfies (see (1.1)- 
(1-2)) 

( _ A + ^ + W^) V = ^ + , (2.5) 
v \x\ ' \x\ 

with WW, W^ 2 \ F«, real analytic in SI C M 3 . 

Assume without loss that Q = Bs(0,r) for some r > 0. (Here, and 
in the sequel, B n (x ,r) = {x E W 1 \ \x — x \ < r}.) Since ip G L 2 (Q), 
Remark 3.2 in Section 3 below implies that ifK is well-defined, as an 
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element of L 2 (K ^\y\ 2 dy). We will show that (p K satisfies (in the 

distributional sense) 

( - A, + 4(W1 1} + \y\ 2 wP))<p K = 4(4 1} + |y| 2 4 2) ) (2.6) 

in K-^Vl) = B 4 (0,y/¥), with Wk, i = 1,2, defined as in (2.3). 
Since W^\F^\ i — 1,2, are real analytic in B 3 (0, r) by assumption, 
and K : R 4 — * M? (see (2.1)) and y \— > ||/| 2 are real analytic, the coef- 
ficients in the elliptic equation in (2.6) are real analytic in -84(0, y/r). 
It follows from elliptic regularity for equations with real analytic co- 
efficients [18, Section 7.5, pp. 177-180] that <p K : 54(0, y/r) -> C is 
real analytic. The statement of Theorem 1.1 then follows from Propo- 
sition 4.1 in Section 4 below. 

It therefore remains to prove that ipx satisfies (2.6). 

By elliptic regularity, y? G W 2 ' 2 (£l') for all fi' = B 3 (0,r'), r' < r. 
(To see this, use Hardy's inequality [28, Lemma p. 169] and that tp G 
W l > 2 (Vt) to conclude that Ay? = G with G G L 2 (fi'). Then use [11, 
Theorem 8.8]). 

It follows that both (Atp) K and (| • \~ 1 tp)x are well-defined, as ele- 
ments of L 2 (.fr _1 (Jl'), ^\y\ 2 dy) (see Remark 3.2 in Section 3 below; see 
also (3.5)). This and (2.5) imply that, as functions in L 2 (i^ _1 (f2')) = 
L 2 (B 4 (0, v^)), 

\y\(Atp) K = \y\((Wtp) K -F K ), (2.7) 

with 

W ( x ) = ^M. + w^ix) , = + F^(x) . (2.8) 

Let / G C^°(-ft'~ 1 (fi)); then there exists r' < r such that supp(/) C 
K- l {VL r ), Q! := B 3 (0,r'); choose {</? n } n6N C C 0O (fi / ) such that <p n -> tp 
and A<p n — >■ Ay? in L 2 (f2')-norm. This is possible since tp G W^ 2,2 (fi'). 
Note that both A/ and 4|y| 2 / belong to C^if -1 ^)) when / does. 
Using (2.4) for y? n G C 0O (Q / ), Remark 3.2 in Section 3 below therefore 
implies that 

/ (Af)(y) tp K (y) dy = Km [ (A/)(y) (<p n ) K (y) dy 
JK-^{n) n ^°°./K-i(n) 

= lim / /(y) [A 2/ ((/? n ) Js: ](y) <iy = lim / A\y\ 2 f(y) (A x tp n ) K (y) dy 
= [ 4\y\ 2 f(y)(A x p) K (y)dy. 
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It follows from this and (2.7) that 

/ (Af)(y)<p K (y)dy = [ A\y\ 2 f{y) {{W V ) K - F K ){y) dy . 

Since (Wip) K = W K ip K , and, by (2.8) and (2.2), 

W K (y) = \y\- 2 (w£\y) + \y\ 2 wP(y)), 
F K (y) = \y\- 2 (F£\y) + \y\ 2 F ( ?\y)), 

this implies that, for all / G C™(K-\Q)), 

VK (y)[-A y f(y) + 4(w£\y) + \y\ 2 wg\y))f(y)} dy 



L 



4(F { K\y) + \y\ 2 Fg\y))f(y)dy, 



which means that (fx satisfies (2.6) in the distributional sense, in 

2.2. The iV-particle problem. In this section we prove Theorem 1.4. 
We only prove the statement (1.12), the proof of (1.13) is completely 
analogous, after an orthogonal transformation of coordinates. We as- 
sume k — 1, the proof for other /c's is the same. 

Let H be given by (1.5). Then, with (x, x') = (x 1} x') G R 3 x R 3JV ~ 3 , 
x' = (x 2 , ■ ■ .,x N ), 

H - E = -A x - A x , - A + V E (x, x') , (2.9) 

\x\ 

where 

N 

v E (x l7 x') = J2-n+ E 1 \- E ( 2 - 10 ) 

3=2 1 ■ 7i l<i<j<N ' 1 31 

is real analytic on Q \ S' x (see (1.8) for S' x ). 
Assume ip G W 1,2 (Vi) satisfies 

(H-E)ip = on Q, (2.11) 

and let (xo,x' ) G fl Ei; then (see (1.10)) xq = 0. We will first 
prove that there exists a neighbourhood Qi(P) of P = (0,x' ) and real 
analytic functions ip^p , i]/p : Qi (P) — > C such that 

^(x) = ^ 1} (x) + |a#j, 2) (x) , x G fix(P) . (2.12) 
By the above, Ve is real analytic on a neighbourhood of (0, x' ), say, 

on 

[/(#) = {(x,x') G R 3 x R 3JV ~ 3 | |x| < R, \x' - x' \ < R} C SI 
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for some R > 0, R small. Let 

U K {R) := {(y,x) 6K 4 x R 3N ~ 3 | \y\ < ^R, \x' - x' \ < R} . (2.13) 

Define now, with K : M 4 -> M. 3 as in (2.1), 

u:U K (R)^C, u(y,x') :=^(K(y),x'), (2.14) 
W : U K (R) -> R , W(y, x') := x') . (2.15) 

Since (by (2.2)) (K(y),x') G U(R) for (y,x') G E/jr(.R), it follows that 
u and are well-defined, and W is real analytic on Uk{R) since K is 
real analytic and Ve is real analytic on U(R). 

As in the proof of Theorem 1.1, we get that (2.11) implies that u 
satisfies 

Q{y,x',D y ,D x ,)u = on U K (R) , (2.16) 

where 

Q(y, x', D y , D x ,) := -A y - 4\y\ 2 A x , + 4\y\ 2 W(y, x') - AZ (2.17) 

is a degenerate elliptic operator, a so-called 'Grusin-type operator'. 

Since \y\ 2 W(y, x') is real analytic on Uk(R), the operator Q is (real) 
analytic hypoelliptic due to [12, Theorem 5.1]. Therefore (2.16) implies 
that u is real analytic on some neighbourhood of (0, x' ) G M 4 x M 3Ar ~ 3 . 

It follows from Proposition 4.4 in Section 4 below that there exist 
a neighbourhood Q^P) C R 3N of P = (0,x' o ) G R 3 x R 3N ~ 3 and real 
analytic functions ipp\ipp^ : ^i(P) — > C such that (2.12) holds. 

Let now 

Q 1 ■= (J Q^P) Gil CM. 3N , 

and define Vi^P = ^i -> C by 

^»(x)=V?(x) when x G fii(P) (i = 1,2). (2.18) 

To see that this is well-defined, we need to verify that if x G Qi{P) H 
fii(Q), then ^ } (x) = ^ } (x), i = 1, 2. Let therefore ^« = - ip%\ 
i — 1, 2, then 

V^ (1) (x) + |a# (2) (x) = , x G fti(P) n fii(Q) , (2.19) 

with i/j^,^ real analytic in Qi(P) n fii(Q). Let x = (0,x{,) G 
f^i(P) C\Qi(Q). Then, since = 1,2, are real analytic, there exist 
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5 > and Pn\i = 1,2, homogeneous polynomials of degree n such 
that 

oo 

$®(x)=^P®(x,J-3> ), * = 1,2, (2.20) 

n=0 

for x G -B 3 jv(x ,5). It follows from (2.19) (by homogeneity) that, for 
all n G N and x = (x, x') G -83^ (x , 5), 

P]l\x,x — x ) + \x\P^2i(x, x — x ) = . 

But for n even, P^ is an even function, while P^-i, an d therefore 
M-f^-i) is odd. Therefore, Prp = P^-i — 0. Similarly for n odd. 
It follows that ip^ = ip^ = on i?3jv(xo,5), and therefore also on 
fti(P) fl Qi(Q), since ip^\ ip^ are real analytic. This proves that ip^ 
and ipf^ in (2.18) are well-defined. Since they are obviously real ana- 
lytic, this finishes the proof of Thereom 1.4. □ 

3. The Kustaanheimo-Stiefel transform 

The KS-transform turns out to be a very useful and natural tool 
for the investigation of Schrodinger equations with Coulombic interac- 
tions. In particular (2.2) and the following lemma are important for 
our proofs. Most of the facts stated here are well-known (see e. g. [13, 
Appendix A]). 

Lemma 3.1. Let K : R 4 -> R 3 be defined as in (2.1), let f : R 3 -> C 
be any C 2 -function, and define fx '■ R 4 — *• C by (2.3). Finally, let 

d d d d 

L(y,D y )-=y 1 - — 2/4^ — h 2/3^ — 2/2^—- (3.1) 

dy A dyi dy 2 dy 3 

(a) Then, with [A; B] = AB — BA the commutator of A and B , 

L(y,D y )f K = 0, [A;L(y,Dy)] — 0, (3.2) 

and (2.4) holds. 

(b) Furthermore, for a function g G C^M 4 ), the following two state- 
ments are equivalent: 

(i) There exists a function f : M 3 — > C such that g = fx- 

(ii) The function g satisfies 

Lg = 0. (3.3) 
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(c) Finally, let U = B 3 {0, r) C M 3 forr G (0, oo]. Then, for (j> G C (M 3 ) 
(continuous with compact support), 

[ \<j ) (K(y))\ 2 dy = U l -^fdx. (3.4) 
JK- l (m 4 Ju Fl 



In particular, 



M<Ml2(^-i({/)) = ^IHIi 2 (t/) • ( 3 - 5 ) 



Remark 3.2. By a density argument, the isometry (3.5) allows to 
extend the composition by K given by (2.3) (the pull-back K* by K) 
to a map 

K* : L\U,dx) ^ L 2 (K-\U)4\y\ 2 dy) 

h-> X 

in the case when [/ = _B 3 (0,r),r G (0, oo]. This makes 0^ well-defined 
for any G L 2 (U). Furthermore, if n — > in L 2 (U), then, for all 
5 G C^X- 1 ^)) (<? G CHX-^C/)), if r = oo) 

lim / g(y)((j) n ) K (y)dy = / g(y)<f) K (y) dy . (3.6) 
^^Jk-uu) Jk-hu) 



>K- l (JJ) JK- 1 ^) 

This follows from Schwarz' inequality and (3.5), 

/ g(y) (Onk(y) - My)) dy 

Jk-i{u) 

<(/ wop dy f mMK -* K 

k Jk--l(u) \y\ J 

f \g(y)\ 2 , y/ 2 , 



LHK-^U)) 



2 



f j M^L^) 1/2 || 0ri _ || L2(l/) ^ O , n^oo. 

^jK-^iu) \y\ 7 



Here the ^/-integral clearly converges since g G C°°(1R 4 ) (g G C^°(1R 4 ), 
if r = oo). 

Remark 3.3. As a consequence of (2.2) and (3.2) (choose f(x) = 
we have that 

L(y,D y )\y\^ = 0, j'gN. (3.7) 

Proof of Lemma 3.1. The lemma is easier to prove in 'double polar 
coordinates' in R 4 . Let 

(R, O) : = (n, r 2 , 9 U 6 2 ) G (0, oo) 2 x [0, 2vr) 2 (3.8) 
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be defined by the relation 

y = y(R, O) = ( yi (R, n),y 2 (R, Q),y 3 (R, ty,y 4 (R, Q)) , (3.9) 

(yi,yd = ri(cos0i,sin0i) , (2/3,2/2) = r 2 (cos6l2,sin6l2) . (3.10) 

Then it follows directly from (2.1) that 



^2 ^,2 

K(,j(RAl))= \ -2r 1 r 2 1 sm(e 1 -e 2 ) | . (3.11) 
2rir 2 cos(6>i - 9 2 ) 

We note in passing that the relation (2.2) is immediate from (3.11). 
In the double polar coordinates, 

' a +£r, (3-12) 



and 



dd 1 d9 2 



, d 2 Id 1 d 2 \ ( d 2 Id Id 2 



sdr\ r\dri r\dQ\) \dr\ r 2 dr 2 r\dQ\ 

Therefore, it is obvious that L and A commute. Furthermore, from 
(3.11) we see that fx only depends on the angles through the expression 
Q\ — 9 2 and therefore, 

The proof of (2.4) is merely an elementary computation, which we leave 
to the reader. This finishes the proof of point (a) of the lemma. 

From (3.2) we infer that in order to prove point (b) we have to show 
that (ii) implies (i). 

To do so, we need to define a function / : IR 3 — > C such that g = fx- 
If x = 0, let f{x) := 0(0), then g(0) = f(0) = f(K(0)) = f K (0) by 
(2.2). Assume now that x G M 3 \ {0}. We claim that the pre-image 
of x under K, K^^x}), is a circle in M 4 (in the literature called the 
'Hopf circle') and that g is constant on this circle. Then, taking any 
y G K~ 1 ({x}) and letting f(x) := g(y), we have that / is well-defined, 
and satisfies fxiv) = f(K(y)) = f(x) — g(y)- This will finish the proof 
of point (b) of the lemma. 

To prove the claim, assume first that x G M 3 \ {0}, x = (xi,x 2 ,xz) 

with (x 2 ,x 3 ) 7^ (0,0). Then the equations (see (3.11) and (2.2)) 

2 2 
r x -r 2 = x u 

—2r\r 2 sin-# = x 2 , 

/ 2 1 2\2 2 1 2 1 2 

(r x + r 2 ) = x l + x 2 + x 3 
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uniquely determine r±, r 2 G (0, oo), and determine $ modulo 2n; choose 
the solution d G [0, 2n). That is, the pre-image of x under K is the 
set of points in M 4 with double polar coordinates (n, r 2 , 0i, 2 ), where 
(ri,r 2 ) is the unique solution to (3.13), and 9\ — 9 2 = i? modulo 2n, 
with i? G [0, 2ir). Defining new angles 6 = 6i + 6 2 ,d = 6\—9 2 , this set is 
the circle in M 4 with centre at the origin and radius (x\ + x 2 + X3) 1 / 4 = 
\Jr\ + r|, parametrized by G [0, 27r). Since, by (3.12), the function 5 
(strictly speaking, g composed with the map in (3.9)) is independent 
of 9 = Q\ + 6 2 , g is, as claimed, constant on this circle. 

On the other hand, assume x = (t,0,0),t G K. \ {0}. Then the 
equations 



(3.14) 



have a unique solution (ri,r 2 ); in fact, (ri,r 2 ) = (y^, 0) if t > and 
(ri, r 2 ) = (0, >/— if t < 0. In both cases, the pre-image of x under K 
is a circle, namely (see also (3.8)) 

C + = {(V* cos 0i, 0,0, Vising)} G M 4 | 9 1 G [0,2tt)} (t > 0) , 
C_ = {(0, sin 2 , cos 2 ,O)} G M 4 | 2 G [0,2tt)} (t > 0) . 

Since y 2 = y 3 = for any y = (yi,y 2 ,ys,y4) G C+, (3.1) and (3.3) 
imply that dg/d0i = 0, with X the angle parametrizing C + , and so g 
is, as claimed, constant on C + ; similarly for C_. This finishes the proof 
of point (b) of the lemma. 

We finish by proving point (c); this is merely a calculation which we 
for simplicity also do in 'double polar coordinates': Recall that \y\ 2 = 
rl + r 2 2 = \x\ (see (2.2)). By (3.10) and (3.11), and since U = B 3 (0,r) 
and K'\U) = £4(0, v^), 



f fin ( fVr r\/r-r\ ^2tt 

/ \<P(K(y))\ 2 dy= \ ndn r 2 dr 2 dQ x 

J K- X {V) JO L JO Jo Jo 

\(f){rj-r 2 2 , -2nr 2 sin(0! - 2 ), 2r x r 2 cos(0 1 - 6 2 )) | 2 } d0 2 . 

In the triple integral inside {•} we make (for fixed 2 ) the change of 
variables 



x = if0 2 (r-i,r 2 ,0i) = [r\-r\, -2r 1 r 2 sin(0 1 - 2 ), 2r x r 2 cos(0i - 2 )) 



ANALYTIC STRUCTURE OF COULOMBIC WAVE FUNCTIONS 



13 



From the foregoing (see after (3.14)) it follows that the image of Kg 2 
is U. The determinant of the Jacobian is 

det(DK e2 ) 

2ri -2r 2 

-2r 2 sin(0i - 2 ) -2r x sin(^ - 2 ) -2rxr 2 cos(0x - 2 ) 

2r 2 cos(^i - 2 ) 2ri cos(0i - 2 ) -2rxr 2 sin(0! - 2 ) 

= 8r 1 r 2 (r2 + r^). 

Recall that |?/| 2 = r\ + r\ — \x\. Therefore the integral is 

/ \ m > dy= f'{! w,)i 2 ^}^ = \ f 

JK-^(U) Jo L Ju °FI J 4 Ju Fl 

This proves (3.4); applying it to a/N"^ gives (3.5). This finishes the 
proof of point (c), and therefore, of Lemma 3.1. □ 

Lemma 3.4. Let the differential operator L = L(y, D) be given by 
(3.1), and let P 2 k be a harmonic, homogeneous polynomial of degree 2k 
in R 4 such that LP 2k = 0. 

Then there exists a harmonic polynomial in IR 3 , Y k , homogeneous of 
degree k, such that 

P 2k (y) = Y k (K(y)) forallyeR\ (3.15) 

with K : R 4 -> R 3 from (2.1). 

Proof. Using that LP 2k = we get from Lemma 3.1 the existence of 
a function Y k such that P 2 k(y) = Y k (K(y)). Since the KS-transform is 
homogeneous of degree 2, Y k is necessarily homogeneous of degree k. 
Furthermore, by (2.4), Y k is harmonic. So we only have left to prove 
that Y k is a polynomial. 

Let C n be the (positive) Laplace-Beltrami operator on the sphere 
gn-i_ Then one can express the Laplace operator in IR n as 

q2 fi — \ jC n 
dr 2 r dr r 2 ^ 
Furthermore, cr(C n ) = {£(£ + n — 1)}^ and the eigenspace correspond- 
ing to the eigenvalue 1(1 + n — 1) is exactly spanned by the restrictions 
to § n_1 of the harmonic, homogeneous polymials in IR ra of degree £. 

Using the fact that AY k = and that Y k is homogeneous of degree 
k in IR 3 we find that Y k \^ 2 is an eigenfunction of £ 3 with eigenvalue 

k(k + 2). Thus there exists a homogeneous, harmonic polynomial Y k 
of degree k such that 
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Since the functions have the same homogeneity, they are identical ev- 
erywhere. This finishes the proof of the lemma. □ 

4. Analyticity and the KS-transform 

In this section we study the regularity of functions given as a com- 
position with the Kustaanheimo-Stiefel transform. We start with the 
one-particle case. 

Proposition 4.1. Let U CK 3 be open with G U, and let ip : U — > C 

be a function. Let U = K~\U) C R 4 , with K : R A -> R 3 from (2.1), 
and suppose that 

ip K = ipoK :U (4.1) 

rea/ analytic. 

Then there exist functions ip^\(p^ 2 \ real analytic in a neighbourhood 
of E R 3 , such that 

(p(x) = V {l \x) + \x\^ 2 \x). (4.2) 

Proof. Note that K(-y) = K(y) for all y G R 4 , so that <p K (-y) = 
(Pk(u) f° r all y G M 4 . It follows that can be written as an ab- 
solutely convergent power series containing only terms of even order. 
Furthermore, since the sum is absolutely convergent, the order of sum- 
mation is unimportant, and so, for some R > 0, eg G C, 

oo 

<p K { y ) = Yl °^ = EE °^ for \y\< R - ( 4 - 3 ) 

/3eN 4 ,|/3|/2eN n=0 \/3\=2n 

This implies (see e. g. [23, sections 2.1-2.2]) that there exists constants 
C 1 ,M 1 > such that 

M < CiMf 1 for all (3 G N 4 . (4.4) 
Note that for fixed n G N, 

Q {2n \y)--= E c ^ ( 4 - 5 ) 

/3eN 4 ,|/3|=2n 

is a homogeneous polynomial of degree In. By [29, Theorem 2.1], 

Q {2n \y) = Y.\y\ 2j H?:l 2] {y), (4.6) 
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where H^} 2 j * s a homogeneous harmonic polynomial of degree 2n — 2j, 
j — 0, 1, . . . , n. It follows that 

oo oo n 

<p K (y) = Q (2n) (y) = E E \y\*H£> v (y) ■ (4.7) 

n=0 n=0 3=0 

We need the following lemma. 

Lemma 4.2. There exist harmonic polynomials Y^f : M 3 — > C, /io- 

mogeneous of degree n — j , such that 

H?:Uy) = Y™(K(y)) for all yeR\ (4.8) 
with K : IR 4 — ► IR 3 /rom (2.1). In particular, the function 

n 

q^Hx):=£\x\*Y™(x) ( 4 - 9 ) 
g (2n) (^(y)) = g (2n) (y) /or a// y G M 4 . (4.10) 



satisfies 



Proof of Lemma \.2 : Recall (see (3.2)) that, with L = L(y,D y ) as in 
(3.1), Ly?x = 0) an d therefore, since power series can be differentiated 
termwise (see (4.7)), 

oo 

= L VK = Y,LQ {2n) ■ (4.H) 

n=0 

Since LQ^ is again a homogeneous polynomial of degree 2n, it follows 
that 

LQ (2n) = 0, n = 0,l,... . (4.12) 
Since L is a first order differential operator, (3.7) implies that 

*i 2 ^_y = \v\ 23 [lh?:i 2j \ , (4.i3) 

where LH^}_ 2 - is again a homogeneous polynomial of order 2(n — j). 
Then (4.6), (4.12), and (4.13) imply that 

n 

E M 2 ^#£ n _y (y) = for all yeR*. (4.14) 

3=0 

Since H^} 2 - is harmonic, and (see (3.2)) [A; L] = 0, we get that 

LH 2 ™}_ 2 - is a homogeneous harmonic polynomial of degree 2(n — j). 
Note that for \y\ = 1, the left side of (4.14) is a linear combination of 
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spherical harmonics of different degrees. From the linear independence 
of such spherical harmonics it follows that 

LH { ^} 2j = for all j = 0, . . . , n, and n e N . (4.15) 

From Lemma 3.4 it follows that there exist harmonic polynomials in 
K 3 , i homogeneous of degree n — j, such that (4.8) holds. 

Now (4.10) follows from this and (2.2). This finishes the proof of the 
lemma. □ 

Lemma 4.2, (4.7), and \K{y) \ = \y\ 2 , imply that 

oo n 

tp K { y ) = tp{K{ y )) = Y.Y.\ K (y)\ jY n-h K (y)) ■ ( 4 - 16 ) 

n=0 j=0 

Formally, we can now finish the proof of Proposition 4.1 by defining 

oo n 

n=0 j=0,j even 
oo n 

^ (2) W:=E E M^O*)- ( 4 - 18 ) 

n=0 j=l,j odd 

However, it is not a priori clear that these sums converge and thus 
define real analytic functions. The remainder of the proof will establish 
the necessary convergence. 

Lemma 4.3. There exists r > such that the two series in (4.17) and 
(4.18) converge for \x\ < r. 

More precisely, there exists a universal constant R > such that 
with C x := Rd, Mi = 2Mf (with C u M x from (4.4) ), 

\Y^]{x)\ < C x M™\x\ n ~i. (4.19) 

Proof. Clearly, the convergence of the series in (4.17) and (4.18) is a 
consequence of (4.19): take r < l/(2M 1 ). Thus we only have to prove 
the estimate (4.19). 

We return to (4.3). For fixed f3, with \(3\ = 2n > we have (again 
using [29, Theorem 2.1]) that, for some df ] E C, 

y p = Y.\v\ 23 df ) P^ 2j {y), (4.20) 

3=0 

where -P2n-2j ^ s a harmonic homogeneous polynomial of degree 2n — 2j, 
which depends on (3, and satisfies ||-P2n-2jlU 2 (s 3 ) = 1- ^ follows from 
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(4.5) and (4.20) that 

Q i2n \y) = E \y\ 2j E c?dfpW 2j (y) . (4.21) 

j=0 \/3\=2n 

Comparing (4.6) with (4.21) we see that 

n 

E \y\ 2j [h^v) - E ^f )p ^y)] = • ( 4 -22) 

3=0 \/3\=2n 

Restricting to \y\ = 1, (4.22) becomes a sum of spherical harmonics 
with different degrees, which are linearly independent, implying that 
(see (4.8)) 

Y<$(K(y)) = H?:l 2] {y) = E cpdfP^iy) . (4.23) 

\(3\=2n 

We are now going to bound the Y^^'s in Since the (restriction 

to S 3 of the) P 2 ( fl 2 /s in (4.20) are orthogonal in L 2 (§ 3 ) (they are ho- 
mogeneous of different degrees), we get (by Parseval's identity), from 
setting \y\ = 1 in (4.20), that 

V|df| 2 = f \yP\ 2 du< [ \y\ 2 Wdw= [ l^ = Vol(§ 3 ), (4.24) 

~^ J& i§3 JS3 

and so the d^'s are bounded, uniformly in j and /3, by Vo^S 3 ) 1 / 2 . 

Due to homogeneity, and using [27, Lemma 8], we get, for any y e 
K 3 \ {0} and j < n, that 

= |y| 2n - 2j l^-(2//bl)| < lyl^'Hf^L^sa) 

Vol(S 3 )V2 Vol(S 3 )V2- ^ Zb) 

Note that (see [29, pp. 138-139]) 

#{aGN*|M=€}=(*+^ 1 ), (4.26) 



and so 



= i(2n + 3)(2n + 2)(2n + 1) < 10n 3 . 
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It follows from (4.23), (4.4), (4.24) (4.25), and (4.27) that (with d and 
Mi the constants in (4.4)) 

\y^(K(y))\ < M \df\ KLM (4-28) 

\/3\=2n 

< 10C in 4 \y\ 2n - 2j M^ n = 10dn 4 \R(y)\ n - j M 2n . 

The desired estimate (4.19) clearly follows, using the surjectivity of R, 
with R := 10 max n n A 2~ n . □ 

Recall that each term {x^Y^j (x) in the definition (4.17) of (f^ is 
a homogeneous polynomial (of degree n) in x, and similarly for 
Therefore, the series (4.17) and (4.18) are convergent power series. This 
implies that (p^\(p^ define real analytic functions on {\x\ < r} (see 
[23, sections 2.1-2.2]). 

Finally, using (4.16), (4.17) and (4.18), 

^{K{y)) + \K{y)\^\K{y)) 

oo n 

= J2lL\ K (y)\ jY n-H K (y)) = <p( K (y))- ( 4 - 29 ) 

n=0 j=0 

This, and the surjectivity of R, imply (4.2) and therefore finishes the 
proof of Proposition 4.1. □ 

For the iV-particle case, we have the following analogous result. 

Proposition 4.4. Let U C IR 3 , U' C JR 3JV ~ 3 be open, with e U, 
x' e U' and let : U x U' -> C be a function. Let U = R' 1 ^) C M 4 , 
with R : IR 4 — > IR 3 /rom (2.1), and suppose that 

u : W x [/' -> C (4.30) 
(w,x) i-> ip(R(y),x') 

is real analytic. 

Then there exist functions ^\^ 2 \ real analytic in a neighbourhood 
W of (0,^) G K 37V ; such that 

ip{x, x) = i> {l) (x, x') + \x\i> {2) (x, x') , (a:, x') E W . (4.31) 

Proof. Define 

<Py(x) ■= ^2^(x,x')\ x , =x , , ^ K (V) ■= <Py(K(y)) . (4.32) 
This is well defined by the assumption on u. 
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Since, as in the proof of Proposition 4.1, u is even with respect 
to y G IR 4 , and the series converges absolutely, we have, for \y\ < 
\/R, \x' — x' \ < R for some R > 0, cg 7 G C, 

u (y, x ') = E °^ yP ( x ' ~ ' 

/3GN 4 ,|/3|/2GN, 7 eN 3JV - 3 

with 

M < C 2 Mf +71 = CWf M^ 1 for all (5 G N 4 , 7 G N 3N ~ 3 , (4.33) 
for some constants C 2 , M 2 > 0. Clearly it follows that 

<PiAv)= E c ^/> ( 4 - 34 ) 

/3eN 4 ,|/3|/2eN 

so that 

«(y.^) = E ( E c Pl yp)(x'-x> y 

7e N 3JV - 3 /3eN 4 ,|/3|/2eN 
= E ^7,-K-(j/) ( X> -^o) 7 - ( 4 - 35 ) 

Moreover, from (4.33) we have that, for all 7 G N 3JV ~ 3 , 

M< C^Mf where C 1 ( 7 ) := C 2 M^ 1 . (4.36) 

In particular, (4.34) and (4.36) show that ip lt K is re al analytic near 
y = 0. Repeating the arguments in the proof of Proposition 4.1 for 
(f JtK for fixed 7 G N 3A ^ 3 , we get that 

00 [n/2] 

^) = EEN 2<y S» ( 4 - 37 ) 

n=0 £=0 

00 [(n-l)/2] 

+ N E E ■^n-(2£+l)( X ) ' 

n=0 £=0 

where F^' 7 : M 3 -> C are harmonic polynomials, homogeneous of de- 
gree n — k, depending on 7 G N 3JV ~ 3 . Therefore, for some 0,3,(7), ^(7) G 
C.aeN 3 , ' 

[n/2] 

Elf«» = EM7K, (4.38) 

£=0 |a|=n 
[(n-l)/2] 

E M 2£ ^m)(*)= E M7)*°, ( 4 -39) 

£=0 \a\=n-l 
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with (see (4.19)), 

| E M7K| < i?C 1 ( 7 )n(2M 2 2 ) n |a:r, (4.40) 

\a\=n 

I E Ul> a \ <RCi{l)n(2M%) n \x\ n . (4.41) 

|a|=n— 1 

Recall that (see (4.26)) 

#{ 7 6N 3 ^ 3 || 7 | = A;}=f 3i ! + A; ; 4 V (4.42) 



37V -4 

By definition, discarding part of the denominator, 

This last product contains (3iV — 4) terms each of which are smaller 
than {3N + k). Thus 



3N + k - 4 
3iV-4 



< (3iV + A;) 37V - 4 <C3A; 3Ar , 



for some C3 (depending on iV) and all k > 1. 

It follows that, for |x| < 1/(4M 2 ), \x' - x' \ < 1/(2M 2 ), 

00 

E 

7£ N3JV-3 n =0 |a|= n 



<^E E E( 2M 2 2 ) n i^r M ^ 

fc=o 7 eN 3JV - 3 ,| 7 |=fc «=o 

OO , Q AT OO 

fc=0 n=0 

and so, with a Q7 := a a (7), 

V> (1 W):= E E a ^(^-^o) 7 ( 4 - 4 3) 

T g N 3JV-3 a6N 3 

defines a real analytic function in a neighborhood of (0, x' ). Similarly, 
with b ai := 6 a (7), 

^ 2 )(x,x'):= E E &^ Q (^-4) 7 (4-44) 

7eN 3JV-3 ag pj3 
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defines a real analytic function in a neighborhood of (0, x' ). From the 
above observations and from (4.35), (4.32), (4.37), (4.38) and (4.39) it 
follows that 

^(K(y),x') + \K(y)\^\K(y),x') 

oo [ra/2] 

= E EEi^)i" y S 7 (^(i/))^-4r 

oo [(n-l)/2] 

+ E \ K M\ E E 

7 gN3iv-3 n=0 €=0 

= E ( Pj,i<(y)(x' -x' o y = u(y,x') =i/j(K(y),x'), 

and so, by the surjectivity of K, 

i/>(x,x?) =if>W(x,x') + \x\^ (2 \x,x') , 

with ip( l \i = 1, 2, real analytic on 

{(x,x) G R 3N | |x| < 1/(4M 2 2 ) , \x' -x' \ < 1/(2M 2 )} . 

This finishes the proof of Proposition 4.4. □ 
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